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We develop a method for constructing exact cosmological solutions of the Einstein equations based 
on representing them as a second-order linear differential equation. In particular, the method allows 
using an arbitrary known solution to construct a more general solution parameterized by a set of 37V 
constants, where N is an arbitrary natural number. The large number of free parameters may prove 
useful for constructing a theoretical model that agrees satisfactorily with the results of astronomical 
observations. Cosmological solutions on the Randall- Sundrum brane have similar properties. We 
show that three-parameter solutions in the general case already exhibit inflationary regimes. In 
. contrast to previously studied two-parameter solutions, these three-parameter solutions can describe 

' an exit from inflation without a fine tuning of the parameters and also several consecutive inflationary 

fSJ I regimes. 
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I. INTRODUCTION 

d 

' A method for constructing and analyzing exact cosmological solutions of the Einstein equations based on repre- 

• , senting them as a second-order linear equation (we call this the linearization method in what follows) was presented 

1^ ■ in [H (other methods for constructing exact solutions in cosmology can be found in 0-0) • Indeed, it is easy to see 

II ' that in the case of the of the flat Friedmann metric, the third power of the scale factor t/j — satisfies the equation: 
O 

a 

where p is the density and p is the pressure of the matter filling the universe. Here and hereafter, we use the system 
of units with SttG/S = c = 1. In the case where a minimally coupled scalar field 4> with the self-interaction potential 
^ V{(p) is dominant and in the presence of a cosmological constant with the density A Eq. ([T]) formally coincides with 
the Schrodinger equation 

, where the potential is U{t) — 9V, and the spectral parameter is A = — 9A. In the quantity V is assumed to 
' be a function of time: V{t) = V{(j){t)), which was called the history of the potential in [l|. Giving an explicit form 
\ U{t) together with the corresponding boundary conditions allows finding the general solution of ([2]). An important 
• • ■ consequence of this investigation is that the regime is independent of or weakly dependent on the type of the potential, 
which is quite significant for the whole theory. Unfortunately, the problem of the end of inflation turns out to be 
substantially more difficult, and solving it in the framework of the approach described above apparently involves 
^ ' additional assumptions (the authors of [i| proposed modifying the potentials to make them depend on the temperature. 
The establishment of a Friedmann regime can then be described as a phase transition in the matter of the early 
universe). The study of Eq. ^ in its application to cosmology was continued in where the Darboux 

transformation was used to construct new exact solutions (a similar technique was used in [lO| to construct exact 
solutions on the brane and on the encompassing space carrying an orbifold structure). 

Here, we present a modification of this method. The Einstein equations imply that the function ^„ = a", where 
n n is an arbitrary (not necessarily integer) number, satisfies a Schrodinger equation with a function J7„, that is a 
linear combination of the density and pressure. If we assume that the universe is filled with a minimally coupled 
scalar field, then [/„ is a linear combination of the potential V{(j}) and the kinetic term, and it is therefore no longer 
reasonable to call C/„ the history of the potential. Hereafter, we call the quantity J7„, the "potential" in quotation 
marks to distinguish it from the self-interaction potential V. 
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Fixing a "potential" we can find solutions for the function i/)„ and thus find the scale factor a„ = -tpl/^ as a function 
of time. In general, a solution of the Schrodinger equation has the form 

M-t)^ciMt) + c2Mt), (3) 

i.e., it depends on two arbitrary constants. The scale factor in turn depends on three parameters: a = a(n, Ci, C2; t). 
This circumstance allows constructing multiparameter solutions of the Einstein equations as follows. We assume that 
the quantity a(t) is determined, for example, by astronomical observation. Raising this function to the nith power, 
we obtain a function denoted by V'l in ([3|). In the next step, we find ^^2 based on the condition of linear independence, 
which we write as 

after which we find a three-parameter solution a = a{ni, c'j^\ c^"^ ;t). We can now repeat this procedure and find a 
six-parameter solution 

t (1) (1) (2) (2) 

and so on. After N steps, we obtain a solution that depends on 3A^ parameters. On the other hand, knowing the 
scale factor in the Friedmann cosmology allows computing all other characteristics, such as the Hubble parameter, the 
acceleration, the density, and the pressure. The photometric distance, which is a crucial quantity for testing models, 
is a function of six parameters; therefore, using the above procedure sufficiently many times, we can make the model 
consistent with observations. 

Unfortunately, the effectiveness of this method is considerably reduced when considering Friedmann models with 
a nonzero spatial curvature. In this case, the function ■(/)„ = a" for an arbitrary n n satisfies not the Schrodinger 
equation but an equation with an additional nonlinear term. There are no meaningful techniques for integrating 
such an equation; it is therefore generally difficult to find a solution depending on two arbitrary constants, which 
is necessary for a complete examination of the problem, similar to the investigation in[l|. It would be interesting 
to show that the conclusions in [l| also hold for fc = ±1. Fortunately, solving the complicated nonlinear equation 
discussed above is unnecessary for this: it suffices to set n — 1 instead of n = 3. In this case, the wave function is the 
scale factor itself. One of the two cosmological Einstein equations in the Friedmann metric is a second-order linear 
equation with a potential proportional (with a minus sign) to p + 3p/c^. The physical meaning of the "potential" is 
also more transparent. If we set the cosmological term to zero, then if the potential is negative, then the strong energy 
condition is satisfied, and if it is positive, then this condition is violated, which generally implies inflation. Therefore, 
the problem of studying inflationary regimes becomes much simpler, even when compared with flat (fc = 0) models. 

This paper is organized as follows. In Sec. 2, we formulate the linearization method exactly, i.e., the reduction of 
the Friedmann equations to a Schrodinger equation for arbitrary n. Moreover, we prove a similar assertion for the 
cosmology on the RandallSundrum I brane (RS-I). In Sec. 3, we give several examples of exact solutions. We see 
that three-parameter families of solutions are much richer in properties than the two-parameter solutions studied in 
In particular, such "potentials" can lead to solutions describing several inflationary stages. Perhaps our universe 
is currently undergoing one of them. 



II. THE LINEARIZATION METHOD 



We consider the Einstein equations in the Friedmann metric: 

— = , (4) 

- - -Up + 3p). (5) 
a 2 

We assume that a = a{t), p = p{t), p = p(t) is a solution of these equations for fc = 0. Then the function tpn — a" is 
a solution of the Schrodinger equation 

i'n = Un[t)lpn, (6) 

where the "potential" is 

C/„(t)=n2p-^(p-(-p). (7) 
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If the universe is filled with a scalar field (f) with the Lagrangian L — ^ — ^(0); then 

U^^'^^^^4>^+n'V{4>). (8) 

Remark 1. For n = 3, the "potential" is — ^V((t){t)). This case was studied in detail in [ij, where this quantity, 
as already noted, was called the history of the potential because appears in the equations as a function of time 
[U = U{t)), not a function of the field variable (j)- Nevertheless, the potential U ^ V, the physical meaning of the 
potential U therefore seems clear. If n ^ 3, then [/„ is a certain linear combination of the kinetic term (j)'^ /2 and the 
self- interaction potential V (. ([7]), dH])), and the physical meaning of such a J7 is not so evident. 

Nevertheless, the effectiveness of the method presented in [l| (and developed in Q, [§]), precisely consists in 
reducing a complex nonlinear problem to a linear equation. This allows nding full two-parameter solutions, which 
exhibit inflationary behavior under very general assumptions. The fact that U basically coincides with V was not used 
anywhere in those papers and therefore played no role. Similarly, we here consider a generalization of this method to 
arbitrary n. 

Furthermore, the physical meaning of the "potential" U3 is clear only for a universe filled with a scalar field. If we 
consider a universe in which, for example, electromagnetic radiation is dominant, the physical meaning of the quantity 
U3 becomes ambiguous. 

Remark 2. If we assume that the universe contains a nonzero vacuum energy with the density p^c'^, in addition 
to the matter fields, then Eq. ^ takes the form of the spectral problem 

i^n = {Unit) - Xn)^n, (9) 

where the spectral parameter is A„ = — 3n^p^/3. Just as for Eq. we can consider a problem for the eigenvalues 
and the eigenfunctions of Eq. ^ if we specify homogeneous initial conditions. As noted in [l[, Eq. has the form 
of a quantum mechanical problem with a discrete spectrum. The fact that each such solution only admits a bounded 
or countable set of allowed values of the cosmological constant (if we specify homogeneous initial conditions) may 
clarify the question of the actual value of the cosmological constant. 

We note that if a solution of Q is known, then we can use ^ and ^ to find the scalar field and the potential: 



We can obtain the dependence V = V{(j)) from these expressions, although it is clearly not always possible to do this 
explicitly. In the general case, a solution of ^ has the form 

^'n = Ci^n + C2^„, (12) 

where ipn is a linearly independent solution with the same potential: 

/■* dt' 

Mt) = Mt) J ^ ^ Mtm. (13) 

Equation (fT2|) allows proving the following assertion. 

Assertion. Let a = a{t) be a solution of ^ and ^ ior k = and the corresponding p and p. Then the three- 
parameter function a„ = a(t; ci, C2, n) of the form 

dt ^ 



a„ - a (^ci + C2 y ^ j , (14) 
is a solution of ([5] and [S]) for a new energy density pn and pressure pn satisfying the condition 

" Pn - Y \Pn + Vnic ) = n p - — [p + p/c ) . (15) 
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Remark 3. In general, this assertion holds for A; = 0. If fc = ±1, then it holds only if n = 0, 1. 
Remark 4. A similar assertion can be formulated for solutions describing the RS-I brane. In this case, the 
Friedmann system is modified by taking the brane tension a into account and becomes 



a J \ 2(7 

a ^^^^ 
-2- = p + 3p + - (2p + ip) , 
a a 



It is easy to see that the function ■(/)„ = a" satisfies the linear Schrodinger equation 



, W^n, (17) 



with the potential 



W„ = ^ (2np-3(p + p) + ^(np-3(p + p)) 



^ . (18) 



2n{K + y) - e/c + -(x + y) {n{K + - m) 

A 



where Y = 1^(0) and K — <fr /2. Hence, the method of generating 3iV-parameter families of solutions described above 
can also be applied to cosmology on a brane. An example of an exact two-parameter solution (with n — 3) was given 
in (ll| , where a linearization method for a simple anisotropic cosmological model was described. 

In what follows, we consider several examples of exact solutions based on an integrable potential of the Schrodinger 
operator, following [l|. 

III. GENERATING EXACT SOLUTIONS WITH A GIVEN Ur,{t) 

We consider the model "potentials" 

Un{t)^tJih^, {A) 

^ cosh^(Aot) ^ ^ 

We study the solutions for potentials (A) and (B) for the possible existence of inflationary regimes and exit from 
inflation. 

Potential (A). For potential (A), the solution of Eq. ^ with zero boundary conditions as t ±oo has the form 

= A"iJ,(/ii) exp(-^iV2), 

where A is a constant and Hg{fit) arc the Hcrmite polynomials of order s. The corresponding evolution of the scale 
factor is 

a{t) = AH^^int) exp{-fit^/2n). (19) 

In the simplest case, we have "ipn = exp(— /ii^/2). In this case, the dependen on n can be eliminated by a simple 
redenition of the parameter ^. This solution was considered in 1], and we therefore do not consider it further. 

Choosing the wave function of the ground state as the solution of problem ^ is not obligatory. We can consider 
solutions in that correspond to excited levels, but these solutions cannot be used on the entire interval on which 
they are defined. The point is that according to the oscillation theorem, the wave function of the sth excited level has 
s zeros, each of which in the cosmological context corresponds to a singularity with the scale factor tending to zero if 
n > 0, and to innity (a Big Rip singularity) if n < 0. We can use solutions for excited levels with numbers s > 1 and 
consider the dynamics described by a part of the eigenfunction on an interval. Such a universe begins and/or ends 
its existence at the corresponding singularity. For example, if we take the function V'n ^ -ffi(pt) exp(— /if^/2) as a 
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solution of ([9]) with the "potential" Un = /i^t^ and let n be -1, then the evolution of the scale factor in this universe 
can be written as 

ait) = ^e^(*^-*o)/2. (20) 

This solution describes a universe in which the scale factor is equal to innity at t = (a Big Rip singularity) , takes its 
minimum value at t = fj,~^^'^ and then begins a never-ending inflationary stage. Using formulas (jlOl andlli p . we can 
also find the asymptotic behavior of the scalar field and the potential for t ~ and as i — s- cx). The kinetic term of 
the energy density of the scalar field is negative, i.e., the evolution described by ([20|) corresponds to a phantom field 
(see, e.g., [TSl-fTs'] for phantom fields). Equation (^0)1 corresponds to a negative cosmological constant A = —3^. 

Potential (B). The solution for potential (B) for A„ ~ — < (which corresponds to a nonnegative value of the 
cosmological constant) has the form 

7/;„ = ci(A - Ao tanh(Aot))e^* + C2(A + Ao tanh(Aot))e-^*. (21) 

If A = Aq, then formula (|2ip simplifies considerably: 



" cosh(Aoi) 

For positive n, the evolution of the scale factor corresponding to the function ipn describes with an inflationary regime 
on the interval {—oo,to), where to is the inflection point of the function cosh^^"(Aot). 

An interesting class of solutions can also be obtained in the simplest case A„ = 0. Then ((2T|) becomes 

■0„ = Ctanh(Aoi). (22) 

If n > 1, then the solution for the scale factor describes a universe leaving the singular state at the moment t = and 
asymptotically approaching a stationary state as t — > oo with a < during the entire evolution. If < n < 1, then 
the universe undergoes an inationary phase until a certain instant and then asymptotically approaches a stationary 
state. 

In the case where A > Aq, ci > and C2 = 0, there is a solution 

ipn = ci(A - Ao tanh(Aot))e'^*, 



a = cy"(A - Ao tanh(Aot))^/"e^*/". (23) 

This solution corresponds to evolution without singularities and is interesting because for certain values of the pa- 
rameters Ao, A and n, formulas (j23p describe a universe that undergoes inflationary expansion on some interval (0, ti), 
then inflation ends, and a secondary inflationary period begins at a time <2 > ti.To determine when such a situation 
occurs, we consider d at the initial instant. Using ([^5)1 . we see that 

a(0) - (1 - n)y'* - 2y2 + 1, 

where we introduce the notation y ~ Xq/X < 1. Therefore, the second derivative of the scale factor is nonnegative at 
the initial instant if 

.,2 ^ „,2 1 



< < % = 



1 



A further investigation shows that satisfaction of this condition means that the universe immediately enters an 
inflationary regime if n > 1. If > ?/g, then solution (|23p describes an evolution with an inflationary regime starting 
at a time io > 0. If rt < 1 then in the narrow interval y^ — A < y^ < yQ, where A ^ yg, there is a period of 
noninflationary expansion, after which the universe again enters a period of accelerating expansion. For example, if 
n = 0.25 and Aq = 0.8A, then the universe undergoes accelerating expansion in the interval (0,0.08/A), followed by 
a noninflationary expansion stage of length approximately 0.4/A. For t > 5/X, the scale factor can be considered to 
change according the exponential law a ~ e"*^^'. 

We conclude by considering the asymptotic behavior of the potential of the scalar field during the early stage of 
the evolution of such a universe. We limit ourself to terms linear in time in the decomposition. Using pH)) and (|lip . 
we find that for i ^ 1/A, 

<Pit) « 00 ± \[^y>^t^ 
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,2 



2(3 -n) 4 

— 2 — y 




Therefore, if there is an initial inflation, then it corresponds to a slow-roll regime. The potential slowly decreases as 
the scalar eld increases or decreases linearly. 



We have proposed a relatively simple method for generating exact solutions of the Einstein- Friedmann equations. 
We used mathematical transformations to reduce the problem to solving the Schrodinger equation for the function 
a" with a "potential" proportional to n^p~ 3n{p + p/c'^)/2. We focused on studying the inflationary regime and the 
exit from it. 

Three-parameter families of solutions exhibit a much richer repertoire of behaviors than do the two-parameter 
solutions studied in flj. Nevertheless, these solutions, like the solutions described in fT], have inflationary phases 
under quite general assumptions. This is an indication that inflation is not something exotic found only in a limited 
number of models. On the contrary, an inflationary regime seems a fairly common occurrence in cosmology not 
requiring any special initial assumptions. 

Unlike two-parameter solutions, three-parameter solutions can have several consecutive inflationary phases, i.e., 
they can not only describe inflation but also describe an exit from the inflationary phase without a special fine tuning 
of the parameters. Moreover, the existence of such solutions can be seen as indirect evidence for the existence of a 
unified realistic model that contains not only an inflationary phase (with an exit from it without a fine tuning) in the 
early universe stage but also a later infiation, which our own universe is perhaps undergoing currently. 

Finally, the procedure for constructing a three-parameter solution can be repeated an arbitrary number TV of times, 
yielding 37V-parameter solutions. The presence of a sufficiently large number of free parameters can be considered a 
defect of a theory but, on the contrary, can be useful for fitting the theoretical model to observational data. We plan 
to return to this question in future investigations. Hence, the described method may prove quite fruitful in cosmology. 
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